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NONLINEAR  EVOLUTION  OP  PLASMA  ENHANCEMENTS 
IN  THE  AURORAL  IONOSPHERE  I: 

LONG  WAVELENGTH  IRREGULARITIES 


1.  INTRODUCTION 

Recently,  large  scale  equatorward  convecting  plasma  enhancements  in  the 
diffuse  auroral  F-region  ionosphere  have  been  identified  and  studied  [Vickrey 
et  al. ,  1980]  using  both  radar  and  satellite  measurements.  Observed  in 
regions  of  diffuse  auroral  particle  precipitation  and  associated  field  aligned 
currents,  these  enhancements  have  overall  latitudinal  dimensions  of  a  few 
hundred  kilometers,  contain  relatively  steep  poleward  and  equatorward  edges, 
and  have  been  shown  to  be  approximately  field-aligned  resembling  vertical 
slabs  of  ionization.  Their  occurrence,  which  is  maximized  in  the  evening- 
midnight  sector,  is  apparently  not  strongly  related  to  magnetic  activity  nor 
to  E-region  processes.  The  presence  of  plasma  density  irregularities  asso¬ 
ciated  with  these  enhancements  has  been  verified  using  satellite  scintilla¬ 
tion  studies  [Fremouw  et  al.,  1977;  Rino  et  al . ,  1978;  Vickrey  et  al.,  1980]. 

The  scintillation  data  have  indicated  that  the  electron  density  irregulari¬ 
ties  are  structured  like  L-shell  aligned  sheets  [Fremouw  et  al.,  1977;  Rino 
et  al.,  1978] .  In  addition,  Rino  and  Matthews  [1979]  have  shown  that  the 
scintillation  enhancements  resulting  from  these  irregularities  cannot  be  ex¬ 
plained  in  terms  of  a  geometrical  enhancement  alone.  A  purely  geometrical 
enhancement  occurs  when  the  signal  propagation  path  intercepts  an  axis  trans¬ 
verse  to  the  magnetic  field  along  which  axis  the  irregularities  have  a  high 

degree  of  spatial  coherence.  Moreover,  the  source  region  of  these  scintilla¬ 
tion  causing  irregularities  has  been  demonstrated  to  be  latitude  limited 

[Rino  and  Owen,  1980]  and  contained  in  a  vertical  slab  of  F  region  plasma. 

Since  these  ionization  enhancements  are  observed  to  convect  equatorward, 

their  poleward  edges  could  be  unstable  to  the  E  x  B  gradient  drift  instabil¬ 
ity  [Simon,  1963;  Linson  and  Workman,  1970]  as  observed  in  artificial  iono¬ 
spheric  plasma  clouds.  Indeed,  for  observed  [Vickrey  et  al.,  1980]  poleward 
Manuscript  submitted  July  16, 1981.  1 
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density  gradient  scale  lengths  of  L  =  10-50  km  and  convection  velocities  of 
approximately  200  m/sec  (Eq  -  10  mV/m)  reasonable  growth  rates  for  the  E  x  IJ 
gradient  drift  instability  can  be  expected  since  y  1=:(BL/cEo)f»  50-250  sec  where 
Y  is  the  E  x  ji  growth  rate,  B  is  the  ambient  magnetic  field  and  c  is  the  speed 
of  light.  Moreover,  it  has  been  shown  [Ossakow  and  Chaturvedi,  1979]  that  by 
applying  the  current  convective  instability  [Lehnert,  1958;  Kadomtsev  and 
Nedospasov,  1960]  the  12  x  B  stable  equatorward  side  of  the  plasma  enhancements 
can  be  driven  unstable  by  the  ambient  field  aligned  particle  precipitation 
currents  in  conjunction  with  the  equatorward  density  gradients.  Other  mecha¬ 
nisms  that  night  account  for  these  irregularities  are  structured  low  energy 
particle  precipitation  and  irregular  field  aligned  currents.  Keskinen  et  al. 
[1980]  showed  that  the  nonlinear  state  of  the  irregularities  in  the  equator- 
ward  edges  of  these  plasma  enhancements  could  be  characterized  by  poleward 
convecting  plasma  depletions  and  equatorward-moving  enhancements.  In  addition, 
it  was  demonstrated  that  these  irregularities  could  be  characterized  by  in¬ 
verse  power  laws  in  the  nonlinear  regime.  However,  these  studies  addressed 
only  the  linear  and  nonlinear  evolution  of  the  equatorward  side  of  the  plasma 
enhancements  and  did  not  include  the  E  x  B  unstable  poleward  edge. 

In  this  report  we  study  the  stability  and  nonlinear  evolution  of  "two 
sided"  models  of  plasma  enhancements  initially  latitudinally  confined  in  order 
to  provide  a  more  realistic  picture  of  the  evolution  of  ionization  enhance¬ 
ments  in  the  auroral  F  region  ionosphere.  In  Section  2  we  present  a  linear 
stability  analysis  of  the  plasma  fluid  equations  which  describe  the  evolution 
of  plasma  enhancements  in  the  auroral  F  region  ionosphere.  The  effects  of  am¬ 
bient  auroral  electric  fields  of  arbitrary  magnitude  and  direction  are  in¬ 
cluded.  In  Section  3  we  describe  the  numerical  methods  used  to  solve  these 
equations.  The  results  of  these  simulations  are  presented  in  Section  4  and  a 
discussion  of  these  results  is  given  in  Section  5. 
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2.  EQUATIONS  OF  MOTION  AND  LINEAR  THEORY 

For  wavelengths  greater  than  the  ion  mean  free  path  we  use  fluid 
equations  t<~  escribe  the  ion  and  electron  plasma.  The  following  geometry 
is  used:  the  y-axis  is  in  the  north-south  direction,  the  x-axis  points 
west,  and  the  z-axis  is  downward  along  the  magnetic  field.  In  this  report 
we  ignore  the  vertical  density  gradient  which  is  weaker  than  the  horizontal 
plasma  density  gradients  [Vickrey  et  al.,  1980]  in  the  typical  diffuse 
auroral  plasma  enhancements.  The  ion  and  electron  fluids  then  obey  the  fol¬ 
lowing  equations  fChaturvedi  and  Ossakow,  1979] : 
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Here  n  (a  ”  i  or  e)  is  the  species  density  and  I£  is  the  total  electric 
field.  Since  we  will  be  interested  in  low  frequencies  and  long  wavelengths, 
we  have  ignored  inertial  terms  in  the  electron  and  ion  momentum  equations 
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(3)  and  (4).  Equation  (5)  results  from  the  assumption  of  quasineutral 

fluctuations  n  =“  n.  s  n.  In  addition,  v  and  V  refer  to  the  electron  and 

ion  velocities  along  the  magnetic  field  giving  rise  to  the  diffuse  auroral 

current.  The  symbol  v^n  denotes  the  ion-neutral  collision  frequency,  the 

electron-collision  frequency,  c  the  speed  of  light,  Te  “  »  T  the  species 

temperature,  c  the  ion  acoustic  speed  and  Cl .  (Q  )  the  ion  (electron)  gyro- 
S  X  6 

frequency.  We  have  neglected  vgn  compared  with  and  taken  va^a  <<;  1  for 
a  =  i,e  (F  region  approximation). 

Any  two  of  equations  (1),  (2),  and  (5)  provide  a  complete  description  of 

the  problem.  We  will  use  the  ion  continuity  equation  (1)  and  (5).  After 

separating  the  total  electric  field  into  an  ambient  and  fluctuating  part 

E  =  -  V^6cp  and  transforming  to  a  frame  drifting  with  velocity 

V  =  -  (c/B)  [z  x  E  -  (v.  /Q.)E  ]  we  can  write 
— o  —  — o  in  l  — o 


x 


y^dcp'V^n  - 


(vin/Qi)VJ.6<p*V]  = 


n 


n.  n  \ 

_ e_  d  n  j 

v,  v  ,  dz  2j 
in  ei  ' 


(6) 


(7) 


where  =  £^v0  -  VQ) •  Linearizing  (6)  and  (7)  by  separating  n  =  nQ(y)  +  6n 
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where  we  have  assumed  k  =  kx  +  ky"£k  cos  a  +  y  k  sina  and  E  =E  cosBx  + 

— Xx—  y  x  x  — o  o  — 
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in  eq.  (8)  is  maximized  for  k-vectors  propagating  at  angles  a  satisfying 
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where  £  =  (k  /k  )  (a . /v.  )(BV,/cE  )  and  we  have  taken  k  /k  to  be  fixed.  For 
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at  a  density  of  nQ  =  105cm_?)and  Eq  =  10  mV/m  (cEq/B  **•  200  m/sec)  we  find 
£  =“  0.3.  To  lowest  order  for  fixed  8  we  find  from  (9)  the  result  that 
a  =“  p/2  +£■  sin(£/2).  In  other  words,  for  P  ^  0,  the  linear  growth  rate  maxi¬ 
mizes  away  from  the  direction  perpendicular  to  the  initial  density  gradient 
by  the  angle  a  ="  P/2  +  ^  sin  (P/2).  As  a  result  the  maximum  growth  rate  from  (8) 
can  be  written  to  lowest  order 
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For  ambient  electric  fields  perpendicular  to  the  initial  density  gradient 
(P  -  0),  eq.  (8)  implies  that  the  maximum  growth  rate  occurs  for  a  »  0,  i.e., 
k  along  E  x.  This  can  be  written 
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where  0  =  k^/k^  Note  that  by  comparing  eq.  (10)  and  (11)  the  linear  growth 


rate,  for  arbitrary  k,  is  reduced  when  is  not  exactly  perpendicular  to 


VnQ,  i.e.,  p  ^  0  because  of  the  cos2  (g/2)  factor  in  (10).  In  regions  of  the 


plasma  enhancements  where  dn  /dy  <  0  and  for  wavelengths  whose  perpendicular 

o 


and  parallel  diffusive  damping  are  negligible  we  find  the  condition  for 
unstable  growth  [  (v^/Q^)  (cEox/B)  +  Slv^l]  >  0  where  we  have  taken,  for 


example,  the  currents  to  be  downward,  i.e.,  <  0.  For  westward  electric 


fields  Eqx  >  (B/c) (Ci/vin)|0|  |v^|,  the  effects  of  the-  field-aligned  currents 
will  be  to  reduce  (0  <  0)  or  enhance  (9  >  0)  the  E  x  B  instability  growth 
rate.  However,  when  dnQ/dy  >  0  the  condition  for  unstable  growth  becomes 
l  (v^/fl^)  (cEqx/B)  +  0 J  Vd |  ]  <  0  and  could  be  satisfied  for  large  enough  cur¬ 

rents  with  |  ATj|  >  (v^/Qp  (cEqx/b|  9|  )  if  0  <  0.  The  expression  for  the 
growth  rate  y  in  equation  (11)  can  be  maximized  as  a  function  of  0  =  ky/k^,  a 
measure  of  field-alignment,  using  dy/SSL _0  =  0  giving 

V 


v.  cE 

0  =_±n_°*± 


m  B  V. 

1  a 


WfrJ  *  (W) 


(12) 


Using  typical  diffuse  auroral  F  region  parameters  v.  /ft. 
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Eqx  -  10  mV/m,  j. 


n^^eVj  1  p,A/m2,  B  **  0.5G,  nQ  “  105cm-3  this  gives 


| ©nil  10  4»  i*e'»  approximate  field  alignment.  Inserting  these  parameters 


into  eq.  (11)  with  L  =**  20  km,  D  =“0.2  m2/sec  and  Dn  10®m2/sec  we  find  that 


the  fastest  growing  linear  modes  have  growth  times  v-1  “  102  sec. 


3.  NUMERICAL  SIMULATIONS 


Equations  (6)  and  (7)  can  be  written  in  dimensionless  form  by 

introducing  the  following  scaled  quantities  n  =  no/NQ,  6cp  =  6cp/BL,  x  =  x/L, 
*** 

y  =  y/L,  z  =  z/L,  t  =  ct/L  as  follows  (where  we  have  dropped  the  tilde  for 
clarity) 
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with  i  =  1,  ....  9  dimensionless  constants  given  by  Cj  = 

c2  .  (*!„«!>  (V.BL)  +  <vel/n.)(cs2/aiCL),  c,  -  c//vlncL,  S  -  DQ,/.,.,, 

Cs  '  Eox/B'  e!'  *  Eoy/a.  c7  -  O^/v^O^/c),  c8  -  T/eBL,  c9  -  (S^/v^)  Cj. 

In  the  following  numerical  simulations  we  take  advantage  of  the  fact 
that  the  fastest  growing,  most  dangerous  modes  from  linear  theory  are  almost 
field-aligned,  i.e,,  k||/k±  <<  1  where  kjj  (k^)  is  the  component  of  k  parallel 
(perpendicular)  to  the  magnetic  field.  These  waves  are  of  most  interest  to 
us  and,  as  a  result,  we  solve  equations  (13)  and  (14)  in  a  plane  containing 
these  modes  which  is  nearly  perpendicular  to  the  magnetic  field  while  fixing 
the  value  of  k||/kx  1  •  A  similar  approach  has  been  adopted  in  numerical 
studies  of  drift-wave  fLee  and  Okuda,  1976]  and  trapped-particle  [Matsuda  and 
Okuda ,  1976]  instabilities  in  laboratory  plasmas.  The  system  of  equations 
(13)  and  (14)  was  first  transformed  to  the  x  y  z  coordinate  system  (as  shown 
in  Fig.  1)  by  a  simple  rotation  about  the  y-axis  by  the  angle  0  =  kjj/k^  «  1 
using 
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where  0  is  the  angle  for  maximum  linear  growth  rate  defined  by  eq.  (12)  for 

a  definite  set  of  a  parameters  v.  /fi. ,  cE  /BV,,  v  . /Q  .  Since  0  «  1  this 

in  l  ox  d  ei  e 

transformation  can  be  written  d/dx  =-  S/Sx'  ,  b/bz  **•  0  S/dx'  ,  d/d y  =  b/by' 
with  b/bz  ^  0.  As  a  consequence  the  three  dimensional  problem  is  reduced  to 
two-dimensions.  By  solving  equations  (13)  and  (14)  in  the  k  y'  z'  coordinate 
system  a  small  but  finite  k jj  is  effectively  introduced  into  the  model. 

Equations  (13)  and  (14)  were  then  solved  numerically  on  a  mesh 
consisting  of  258  grid  points  in  the  north-south  direction  (y-direction)  and 
102  grid  points  in  the  east-west  direction  (x-direction)  with  constant  grid 
spacing  of  1  km.  As  a  result,  the  simulation  plane,  which  is  taken  to  be 
essentially  horizontal  at  an  altitude  of  350  km  in  the  diffuse  auroral  F 
region,  has  a  north-south  and  east-west  extent  of  256  and  100  km,  respec¬ 
tively.  The  field  aligned  currents  are  taken  to  be  constant  in  space  and 
time  over  the  grid.  The  plasma  density  n  in  equation  (13)  was  advanced  in 
time  using  a  multi-dimensional  flux-corrected  variable  timestep  leapfrog- 
trapezoid  scheme  [Zalesak,  1979]  which  is  second  order  in  time  and  fourth 
order  in  space.  At  each  timestep  the  self-consistent  electrostatic  poten¬ 
tial  of  the  plasma  enhancement  in  eq.  (14)  was  determined  using  a 
Chebychev  iterative  method  [McDonald,  1980]  with  a  convergence  criterion  of 
10"4.  Periodic  boundary  conditions  were  imposed  in  the  east-west  direction 
with  Neumann  boundary  conditions  (3/9 y  =  0)  in  the  north-south  direction.  A  slab 
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approximation  is  used  to  model  the  zero  order  convecting  plasma  enhancements 
in  the  diffuse  auroral  ionosphere  with  the  north-south  profile  given  by 
nQ(y')  =  Nq  {1  +  4.5  [tanh(y'  -yi)/L+  tanh(y2-y')/L)  ]  }  (1  +  e(x',y'))  with 
=  1  x  105cm-3,  yi  =  50  km,  and  y2  =  125  km.  This  gives  a  maximum  plasma 
enhancement  density  to  background  ratio  of  approximately  10.  The  initial  per¬ 
turbation  e(x‘  .y')  has  a  root  mean  square  value  of  0.01  and  has  a  radially 
Gaussian  dependence  on  xy  and  y* .  We  now  drop  the  prime  notation  for  clarity. 

4 .  RESULTS 

In  the  following  we  consider  the  linear  and  nonlinear  evolution  of 
plasma  enhancements  in  the  diffuse  auroral  F  region  ionosphere  in  an  approxi¬ 
mately  horizontal  plane  at  350  km  altitude  almost  perpendicular  to  the  magne¬ 
tic  field.  We  take  the  following  typical  parameters  [Vickrey  et  al.,  1980; 
Schunk  and  Walker,  1973;  Banks  and  Kockarts,  19731  L  =  20  km,  v.  /ft,  =  2  x 

10~4,  v  ./ft  =  2  x  10-4,  E  =10  mV/m  T  =  T.  =  1000°K  and  J..  =  1  pA/m? 
ei  e  ’ox  e  l  II 

(which  gives  a  current  velocity  of  V^  “*  60  m/sec  with  Nq  °*  1  x  105cm-3) .  In 
addition,  we  assume  that  the  diffuse  auroral  particle  precipitation  current 
Jjj  is  downward  (V^  <  0)  and  spatially  and  temporally  uniform  over  the  entire 

plasma  enhancement.  In  order  to  find  the  location  and  magnitude  of  the  maxi¬ 
mum  linear  growth  rates  to  be  expected  with  this  set  of  parameters  we  first 
compute  6^  =  kji/k^  as  given  in  eq.  (12)  with  V^  =  -|v^|  =  -  60  m/sec.  This 

gives  two  values  for  6  which  are  0+  =  1.4  x  10-3  and  9“  =  -  6.5  x  10-4 .  The 

m 

first-  value  0+  gives  a  maximum  linear  growth  rate  Ymax  “  1*1  x  10  2sec-1  on 

the  poleward  side  (dnQ/dy  <  0)  with  linearly  damped  perturbations  Ymax  ““  3  x 

10~3sec-1  on  the  equatorward  side  (dn  /dy  >  0).  The  second  value  9“  gives 

o 

only  a  marginally  unstable  growth  rate  of  Ymax  ““2.1  x  10  4sec_1  on  the 
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equatorward  side  with  damped  fluctuations  -  2.3  x  10-<+sec-1on  the 

poleward  side.  These  results  agree  with  the  experimental  observations 
fVickrey  et  al.,  1980]  that  the  largest  linear  growth  rates  occur  on  the 
poleward  side  of  the  convecting  plasma  enhancements.  In  this  case  the 
effect  of  the  field-aligned  currents  is  to  enhance  the  E  x  B  gradient-drift 
instability  growth  rate  on  the  poleward  side.  The  currents  are  too  weak  for 
the  cases  studied  observationally  to  give  appreciable  growth  on  the  equator- 
ward  side  of  the  plasma  enhancements.  We  will  then  consider  the  evolution 
of  modes  satisfying  8+  =  kjj/k^  =  1.4  x  10-5. 

We  consider  two  models  with  different  initial  electric  field 
configurations.  Model  1  has  Eqx  =  10  mV/m,  Eq^  =  0  while  Model  2  takes 
Eqx  =  8  mV/m,  E  =  2  raV/m.  Figure  2a-2d  gives  the  evolution  of  the 
plasma  enhancement  using  Model  1  (no  northward  electric  field).  Figure  2a 
shows  the  initial  configuration  which  includes  the  small  perturbation. 

Figure  2b  illustrates  the  linear  regime  of  the  simulations  and  shows  unstable 
growth  on  the  poleward  side  of  the  plasma  enhancement  as  predicted  by  the 
linear  result  given  by  eq.  (li).  One  can  note  the  depletion  jetting  to  the 
equatorward  side  of  the  enhancement  in  analogy  to  the  initial  evolution  of 
the  E  x  IJ  gradient  drift  instability  in  artificial  ionospheric  plasma  clouds 
[Zabusky  et  al. ,  1973;  Scannapieco  et  al.,  1976].  Figure  2c  gives  the 
structure  of  the  plasma  enhancement  at  t  =  1000  sec  and  shows  steepened  fin¬ 
gers  which  are  beginning  to  elongate.  Finally  Figure  2d  displays  the  plasma 
enhancement  at  t  =  1600  sec  in  the  well-developed  nonlinear  regime.  The 
poleward  edges  of  the  principal  fingers  (striations)  have  steepened,  become 
quasi-one  dimensional  and  bifurcated.  The  length  scales  on  Figure  2a-d  are 
distorted  with  the  depletions  longer  and  narrower  than  is  depicted. 


Figure  3a-b  give  sample  one-dimensional  spatial  power  spectra  at 

t  =  1600  sec  both  in  the  east-west  (P(k  ))  and  north-south  (P(k  ))  directions 

x  y 

respectively  for  Model  1.  These  power  spectra  are  defined  as  follows 

P(k  )  =  /dk  P(k  ,k  ) 
x  J  y  x  y 

and 


P(k  )  =  f ik  P(k  ,k  ) 
y  J  x  x  y' 


—  —  1  2 

where  P(k  ,k  )  =  (L  L  )  [6n(k  ,k  )/N  ]  is  the  spectral  density,  6n  =  N-N 

xy  xy  xyo  o 

with  the  peak  plasma  enhancement  density,  and  is  the  area  of  the 

numerical  simulation  plane.  For  both  cases  these  power  spectra  are  well- 
fitted  with  an  inverse  power  law  with  spectral  index  n^  2-2.5  for  2-ir/k^ 
between  approximately  100  and  3  km  and  n^  2  for  2ir/k^  between  256  and  3  km. 

Figure  4a-4d  illustrate  the  evolution  of  the  plasma  enhancement  using 
Model  2  with  a  westward  electric  field  E  =  8  mV/m  together  with  a  north¬ 
ward  component  EQy  =  2  mV/m  using  the  same  initial  conditions  as  in  the 
previous  case  (E  =  0)  at  t  =  0,  550,  1000,  and  1600  sec.  As  in  model  1, 

these  electric  field  components  give  two  values  for  0  which  are 

m 

+  -5  -  -4  + 

9  =  1.7  x  10  and  0  =  -  5.03  x  10  .  Since  0  gives  the  largest  linear 


growth  rate  of  y 


max 


-3  -1 

8.0  x  10  sec  on  the  poleward  side,  we  consider 

•s“5 


in  Model  2  those  modes  satisfying  6  =  kjjk^  =  1.7  x  10  5 .  Figure  4a 

gives  the  initial  configuration  which  is  identical  to  Fig.  2a.  Figure  4b 
shows  the  isodensity  contours  of  the  plasmas  enhancement  at  t  =  550  sec  where 
one  can  note  a  westward  tilt  to  the  fingers  on  the  unstable  poleward  side  and 
decreased  depletion  jetting  to  the  equatorward  side  in  comparison  to  Fig.  2b. 
This  tilt  is  reminiscent  of  ionospheric  plasma  cloud  structuring  [  Perkins 
and  Doles,  1975]  in  ambient  electric  fields  which  are  not  initially  perpen¬ 
dicular  to  the  initial  plasma  cloud  density  gradient.  The  tilt  can  be  ex¬ 
plained,  in  part,  by  referring  to  the  discussion  following  eq.  (9)  which 


11 


states  that  the  linear  growth  rate  maximizes  away  from  the  direction  perpen¬ 
dicular  to  the  initial  plasma  enhancement  density  gradient  when 
E  •  VnQ  f  0.  The  decreased  depletion  jetting  to  the  equatorward  side  of 
the  plasma  enhancement  in  Fig.  4b  as  opposed  to  Fig.  2b  can  be  resolved  by 
noting  than  the  linear  growth  rate  in  Model  2  (3  =“  14°)  is  reduced  from  Model  1 
(3  =  0)  due  to  the  cos2  (3/2)  factor  in  eq.  (10).  Figure  4  c  gives  the 
structure  of  the  plasma  enhancement  at  t  =  1000  sec  for  Model  2  which  also 
appears  slightly  less  developed  than  in  Fig.  2c  which  contains  no  northward 
electric  field  component.  Finally,  Fig.  4d  details  the  plasma  enhancement  at 
t  =  1600  sec  for  Model  2.  One  notes  that  the  eastward  edge  of  the  large 
finger  on  the  westward  side  of  the  grid  is  slightly  stepper  than  the  westward 
edge  (E^  •  Vn^  ^  0) .  This  may  lead  to  L-shell  aligned  east-west  kilometer 
size  structures  due  to  secondary  E  x  B  instabilities  although  the  present 
simulations  do  not  have  adequate  spatial  resolution  to  develop  such  an 
hypothesis.  In  addition,  there  is  more  pronounced  bending  of  the  fingers  in 
Fig.  4d  in  comparison  to  Fig.  2d  and  less  bifurcation  on  the  poleward  tips 
of  the  striations.  As  mentioned  previously,  these  features  can  be  explained, 
in  part,  by  the  small  northward  electric  field  components. 

Figure  5a-b  give  the  one-dimensional  east-west  P(k^)  and  north-south 
P(ky)  spatial  power  spectra  at  t  =  1600  sec  for  Model  2.  The  power  laws  and 
spectral  indices  are  similar  to  Model  1. 

5.  SUMMARY  AND  DISCUSSION 

We  have  studied,  through  numerical  simulations,  the  nonlinear  evolution 
of  plasma  enhancements  in  the  diffuse  auroral  F  region  ionosphere.  We  have 


I 

shown  that  equatorward  convecting  plasma  slabs  initially  limited  in  latitudi¬ 
nal  extent  can  be  destabilized  on  the  poleward  sides  by  a  combination  of  the 
effects  of  convection  and  field  aligned  currents.  These  simulations  indicate 
that  this  destabilization  leads  to  striation-like  structures  (elongated  in 


the  north-south  direction)  which  can  form  and  cascade  from  long  wavelengths 

(~  100  km)  to  shorter  scale  sizes  (~  1  km)  on  the  order  of  an  hour.  The 
one-dimensional  irregularity  spatial  power  spectra  in  the  east-west  direction 

P(k  )  «  k  x,  n  =“  2-2.5,  for  2ir/k  between  100  km  and  3  km  while  in  the 

XXX  X 

-ny 

north-south  direction  P(k  )  cc  k  ,  n  ^  2,  for  2if/k  between  256  km  and 

y  y  y  y 

3  km. 


In  this  paper  we  have  studied  the  quasi-two  dimensional  linear  and 
nonlinear  evolution  of  models  of  plasma  enhancements  in  the  diffuse  auroral 
F  region  ionosphere.  This  has  been  accomplished  by  solving  the  plasma  fluid 
equations  in  a  horizontal  plane  approximately  perpendicular  to  the  magnetic 
field.  The  observed  plasma  enhancements  are  three  dimensional  [Vickrey  et 
al. ,  1980] .  However  the  horizontal  gradients  are  much  steeper 
than  the  vertical  denstiy  gradients  allowing  one  to  approximately 
model  the  plasma  enhancements  by  vertical  slabs.  In  addition,  we  have  not 
included  a  full  spectrum  of  finite  kjj  modes  in  these  simulations.  However, 
since  the  modes  with  maximum  linear  growth  rate  have  k||/k^  «  1 ,  the  import¬ 
ant  structuring  processes  will  occur  in  the  plane  nearly  perpendicular  to  the 
magnetic  field. 

Finally,  we  note  that  we  have  not  addressed  the  source  mechanism  of  the 
plasma  enhancements,  their  coupling  to  other  levels,  e.g..  E-region,  the 
nature  of  the  intermediate  wavelength  irregularities  (X  ~  1  km)  in  the  plasma 
enhancements,  and  the  role  of  neutral  winds.  These  topics  will  be  discussed 


iu  future  studies. 
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yjy^north) 


simulations.  The  x'y'  is  the  simulation  plane. 
;',z  axes  are  coplanar. 
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(c)  (d) 

Fig.  2  —  Real  space  isodensity  contour  plots  of  n(x',y')/N0  for  model  1  at 
(a)  t  =  0  sec,  (b)  t  =  550  sec,  (c)  t  =  1000  sec,  (d)  t  =  1600  sec.  The  y-axis  is 
compressed  by  a  factor  of  2.58.  The  distance  between  tic  marks  in  the  x- 
direction  (y-direction)  is  5  km  (12.8  km).  Eight  contours  are  plotted  in  equal 
increments  of  1.25  beginning  at  1.25.  The  observer  is  looking  upward  along  the 
magnetic  field  lines. 
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Fig.  3  —  One  dimensional  (a)  x  power  spectra  P(kx)  and  (b)  y  power  spectra 
P(ky)  at  t  =  1600  sec  for  model  1.  In  (a)  kPx  =  2jt/100  km'1  while  in  (b)  kjy 
-  2rr  /256  km'1.  The  dots  represent  the  numerical  simulation  results;  the  solid 
curve  is  a  least  squares  fit  to  modes  2-30  in  the  x-direction  and  to  modes  2-80  in 
the  y-direction.  The  units  of  P(kx),  P(ky)  are  km. 
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CIELD  COMMMC 
defense  nuclear  agenct 
KIRTLANO  AFB.  Ml  87115 

oia  attn  fcpr 

director 

Interservice  nuclear  weapons  School 
XIRTLAA®  aFB,  Ml  87115 

Ola  attn  Document  Control 


WInMCCS  System  Engiicering  org 
WASHINGTON,  U.C.  20305 
ola  Attn  k.  Crawford 

LOMMAfC£R/Dl  RECTOR 
ATMOSPHERIC  SCIENCES  LABORATORY. 

U.S.  ARMY  tLECTRONICS  COMMA® 

WHITE  8AICS  MISSILE  KANGE,  NM  88002 
Ola  ATTN  UtLAStO  F.  NILES 


J I RECTOR 

iMU  ADVANCED  lECH  CTR 
UNTSVILLE  OFFICE 
\  0.  BOX  1500 
(UNTSVILLE,  AL  35807 
Ola  Attn  ATC-T 
OIcy  Attn  aTC-0 
Ola  attn  ATC-R 


Melvin  T.  Capps 
w.  Davies 
don  Russ 


program  Manager 
jlD  Program  Office 
5X51  CISENHOWER  AVEfWE 
V.EXAICRIA,  Va  22333 

Ola  attn  DACS-bMT  J. 


Shea 


Chief  C-E  Services  Division 
j.S.  army  Communications  CMD 
Pentagon  Rm  1B269 
WASHINGTON,  U.C.  20310 

OIcy  attn  'Services  Division 
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COMMAfOER 

FRADCOM  Tecwiical  Support  activity 

DEPARTMENT  OF  THE  4RKY 
Fort  Mottoum.  ft.  j.  07703 

ula  httn  DRSEL-NL-rtD  h.  eennet 
OIcy  Attn  DRSEL-PL-ENV  n.  Bofke 
01a  httn  j.  E.  Quigley 

COMMAfOER 

HARRY  DlAMOfO  LABORATORIES 

Department  of  the  Army 
2800  Ponder  Mill  Road 

Adelphi.  ft  20783  m 

CCf*CI-lNNER  Envelope:  Attn:  OtLHDKBH) 
01a  Attn  DELHD-TI  m.  Weiner 
01a  Attn  DELHD-RB  R.  Williams 
Ola  ATTN  DELHHtP  F.  WlMENITZ 
Ola  Attn  DEOC-NP  C.  Moazed 

COMMAfCER 

U.S.  ARMY  COMM-tLEC  CNGRG  [NSTAL  AGY 
Ft.  Huachuca.  A z  85613 

01a  Attn  CCC-£J€0  George  Lane 

CQMMAM3ER 

U.S.  Army  Foreign  Science  &  Tech  Ctr 
220  7th  Street.  NE 
Charlottesville,  Va  22901 
01a  httnDRXST-SD 
01a  httn  r.  JONES 

Commander 

U.S.  Army  Materiel  Dev  &  Readiness  CMD 
5001  Eisenhower  Avenue 
Alexandria,  Va  22333 

01a  httn  DRCLDC  J.  a.  Sender 

COWWCER 

U.S.  army  nuclear  ah>  Chemical  AGENa 
7500  Backlicx  Road 
Bldg  2073 

Springfield,  va  22150 
Ola  Attn  library 

Director 

U.S.  Army  Ballistic  Research  labs 
Aberdeen  Proving  Grouno.  ft  21005 
01a  Attn  Tech  lib  Edward  SAia 

COMhArOER 

U.S.  army  SATCQM  AGENa 
Ft.  Monmouth,  nj  07703 

01a  attn  document  Control 

COMMAfOER 

u.S.  army  Missile  Intelligence  AGENa 
RED5TOY*  Arsenal-  al  35809 
Ola  httn  Jim  Gamble 

Director 

U.S.  army  TRADOC  Systems  Analysis  Activity 
white  SAfts  missile  range,  nm  s80C2 
Ola  httn  hTAA-sa 
Ola  httn  TGC/F .  Payan  Jr  . 

01a  httn  hTAa-tac  lTl  j.  hesse 

COMPWfOER 

naval  electronic  Systems  Cohwc 
WASHINGTON.  D.C.  20360 

01a  httn  NAVALEX  034  T.  hughes 
01a  Attn  Pft  117 
Ola  attn  PME  117— T 
01a  attn  Code  5011 

Cowwoing  Officer 

haval  Intelligence  support  ctr 

4j01  SUITLANO  ROAD.  BLDG.  5 
WASHINGTON,  U.C.  20390 

oia  attn  m.  uubbin  sTic  12 

Ola  attn  n  ISC -50 

oia  attn  Code  brum  j.  Galet 


COMMAfCER 

naval  Ocean  systems  Center 
SAN  DIEGO,  CA  92152 

03a  attn  Code  532  w.  moler 
01a  attn  Code  Q23Q  C.  Baggett 
01a  attn  Code  81  «.  Eastman 

director 

naval  Research  Laboratory 
Washington.  D.C.  20375  . 

01a  Attn  cooe  4700  I.  p.  Coffey  (26  as  if  un.  1  a  if  class) 
01a  attn  Code  4701  jack  D.  brown 

01a  attn  Cooe  4780  Branch  head  (150  as  if  un,  i  cy  if  class) 

01a  Attn  Code  7500 

Ola  httn  Code  7550 

OIcy  attn  Cooe  7580 

01a  Attn  Code  7551 

01a  Attn  Code  7555 

01a  attn  Code  4730  E.  McLean 

OIcy  attn  Code  4187 

COMMAfOER 

Naval  Sea  Systems  Commajc 
Washington,  D.C.  20362 

01a  Attn  Capt  r.  Pitkin 

naval  Space  Surveillance  System 
jah.gr en,  Va  22448 

01a  Attn  Capt  j.  h.  Burton 

Officer-In-Charge 
Naval  Surface  weapons  Center 
White  Oak.  Silver  Spring,  md  20910 
01a  Attn  Code  F31 

Director 

Strategic  Systems  Project  Office 
Department  of  the  Navy 
Washington.  D.C.  20376 
Ola  ATTN  NSP-2141 
01a  attn  NSSP-2722  Fred  wimberly 

COMMAfOER 

Naval  Surface  weapons  Center 
Dahlgren  Laboratory 
DAH.GREN.  VA  22448 

01a  attn  Code  DF-14  r.  Butler 

Office  of  Naval  research 
Arlington,  va  22217 
01a  attn  Code  465 
Ola  attn  Code  461 
01a  Attn  Code  402 
01a  attn  Code  420 
oia  attn  Code  421 

COMMAf€)£R 

aerospace  Defense  Command/DC 
Department  of  the  air  Force 
ENT  AFB.  CO  80912 

Ola  Attn  X  d».  long 

COMMAfOER 

aerospace  Defense  Comfwo/xpd 
Department  of  the  air  Force 
ENT  AFB.  CO  80912 
oia  attnxpdqq 
01a  Attn  XP 

air  Force  Geophysics  laboratory 
HANSCOM  AFB.  MA  01731 

01a  ATTN  OPR  HAROLD  GARDNER 
01a  Arm  OPR-1  James  C.  Ulwick 
01a  attn  LKB  kewcth  S.  w.  Chaff  ion 
01a  Arm  OPR  Alva  T.  Stair 
01a  Arm  PK>  jules  Aarons 
Ola  Arm  PHD  jurgen  buchaj 
Ola  Arm  PHD  John  p.  Mullen 


ff  weapons  Laboratory 
KIRTLANd  AFB.  NM  87117 
OIcy  Attn  SIX 

Ola  attn  CA  Arthjr  h.  Guenther 
OICY  ATTN  NTYC  1LT  G.  KRAJCI 


aftac 

Patricx  AFB,  Pl  i2925 
Ola  ATTN  TF/MAJ  wIL£Y 

oia  attn  tn 

air  Force  avionics  Laboratory 
wright-Patterson  aFB,  oh  45453 
Ola  Attn  AAD  Wade  hunt 
Ola  Attn  AAD  Au_en  Johnson 


Department  of  cnergy 
Albuquerque  Operations  office 
P.  0.  Box  5400 
ALBUQUERQUE,  87115 

Ola  attn  DOC  CON  for  d.  Sherwood 

Department  of  Cnergy 
LIBRARY  ROOM  6-042 
Washington.  D.C.  20645 

Ola  attn  UOC  Con  for  a,  labowitz 

EGSG,  Inc. 
los  Alamos  Division 
P.  0.  BOX  809 
Los  Alamos,  nm  85544 

Ola  attn  DOC  CON  for  j.  Breedlove 


Deputy  Chief  of  Staff 

KESEARCH,  DEVELOPMENT,  i  ACQ 
Department  of  the  Air  Force 
Washington.  D.C.  20530 
Ola  Attn  AFRDC 

headquarters 

Electronic  Systems  Division^ 
Department  of  the  Air  Force 
Hanscom  AFB.  rv\  01731 
Ola  Attn  XR  J.  LEAS 

HEADOUARTERS 

Electronic  Systems  Divisiqn/YSEa 

DEPARTMENT  of  THE  AlR  FORCE 
HANSCOM  AFB,  MA  01732 
Ola  attn  YSEA 

headquarters 

Electronic  Systems  jivision/DC 
Department  of  the  hi r  Force 
hanscom  AFb,  nA  01731 

OIcy  Attn  bCKC.  j.l.  Clark 


university  of  California 

LAWENCE  LIVERMORE  LABORATORY 
P.  0.  Box  808 
LIVERMORE.  CA  94550 


UICY  MT I N  1AA.  UJH  FOR  IECH  INFO  DEPT 

Ola  Attn  DOC  CON  FOR  L-389  R.  Orr 
Ola  ATTN  DOC  CON  FOR  L-31  R.  HAGER 
oia  ATTN  DOC  CON  FOR  L-46  F.  Seward 


LOS  Alamos  Scientific  Laboratory 

P.  0.  Box  1683 

Los  Alamos,  87545 

Ola  Attn  Doc  Con  for  J.  volcott 
Ola  attn  Doc  Con  for  R.  F.  Taschek 
Ola  attn  Doc  Con  for  E.  jofts 
Ola  attn  Doc  Con  for  j.  Malik 
oia  attn  Doc  Con  for  k.  Jeffries 
Ola  attn  Doc  Con  for  J.  Zinn 
Ola  attn  Doc  Con  for  P.  keaton 
Ola  attn  Doc  Con  for  d.  westervel” 

Sandia  laboratories 
P.  0.  Box  5800 
Albuquerque,  NM  87115 


COMMAJCER 

Foreign  Technology  Division.  aFSC 
WRIGHT-PATTERSON  AFb.  UH  45435 
ula  Attn  nICD  library 
Ula  ATTN  ETUH  b.  BALLARD 

COMMAfCER 

home  air  development  Center.  aF SC 
URIFFISS  AFB.  NY  13441 

ula  attn  DOC  Ubrary/Tslu 
01a  Attn  uc£  v.  Coyne 

post  office  Box  92960 
vorldway  Postal  Center 
Los  Angeles.  CA  90009 
(Space  Defense  Systems) 

01a  Attn  SZJ 


Strategic  Air  CommandAFFS 
Offutt  AFB.  nb  o8113 

01a  attn  XPFS  i'Iaj  B,  Stephan 
OIcy  attn  ADWATE  Maj  Bruce  Bauer 
01a  attn  NRT 

01a  attn  DGK  Chief  Scientist 


?.  0.  Box  92960 
worloway  Postal  Center 
los  Angeles,  Ca  30009 

01a  Attn  SKA  (Space  Comm  Systems)  H,  Clay  in 
SAMSO/MN 

Horton  HFB.  Ca  92409 
(Minuteman) 

01a  Attn  MNH.  ltc  Kennedy 
COEfWQCR 

Some  Air  Development  Center.  AF3C 
hanscom  aFB,  Ma  01731 

01a  attn  £EP  a.  lorentzen 


Oia 

Oia 

Ola 

Ola 

Ola 

Ola 


attn  Doc  Con  for  w.  crown 
Attn  Doc  Con  for  A,  Thorwrough 
attn  Doc  Con  for  T,  wright 
attn  Doc  Con  for  D.  Dahlgren 
attn  Doc  Con  for  3141 
attn  Doc  Con  for  Space  Project  uiv 


Sandia  laboratories 
Livermore  laboratory 
P.  0,  aox  %9 
LIVERMORE.  CA  94550 

01a  attn  Doc  Con  for  b.  Itjrpky 
ula  attn  doc  Con  for  t.  cook 


Office  of  Military  Application 
Department  of  Energy 
WASHINGTON,  D.C.  20545 

01a  attn  doc  Con  for  Dr .  Yo  Song 


Other  Government 


Central  Intelligence  Agency 
attn  RD/SI,  Km  5G48,  HQ  bldg 
WASHINGTON,  D.C.  20505 

01a  Attn  0SI/PSID  Rm  5F  19 


Department  of  Commerce 
National  Bureau  of  St  awards 
WASHINGTON,  D.C.  '20234 

(All  Corres:  attn  Sec  Officer  for) 
01a  attn  R.  Moore 


Institute  for  Telecom  sciences 
national  Telecommunications  &  Info  Aomin 
boulder,  Co  80303 

01a  ATTN  A.  JEAN  (UNCLASS  ONLY) 
Ola  ATTN  W.  JTLAUT 
Ola  ATTN  D.  CROMBIE 
Ula  ATTN  L.  BERRY 
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NATIONAL  UCEANIC  i  ATMOSPHERIC  ADMIN 
dWIROWIENTAL  RESEARCH  LABORATORIES 
uEPARTMENT  Of  COMMERCE 

BOULDER'  Co  803u2 

OIcy  attn  R.  Grubb 

Ola  ATTN  AERONCMY  lAB  G.  REID 

Department  of  Defense  Contractors 

aerospace  Corporation 

P.  0.  Box  92957  _ 

Los  Angeles.  Ca  9OOC0 

01a  Attn  1 .  Garpunkel 
Ola  attn  T.  Salmi 
Ola  attn  v.  JOSEPHSON 
Ola  attn  S.  Boner 
Ola  attn  n.  Stoooaeu. 

01a  Attn  D.  Olsen 


analytical  Systems  Engikering  Corp 
5  Old  Concord  Road 

BURLINGTON'  MA  01803 

Ola  attn  Radio  Sciences 

BERKELEY  RESEARCH  ASSOCIATES.  INC. 
P.  0.  BOX  983 
BERKELEY.  CA  94701 

Ola  attn  j.  wrkman 

Boeing  Company,  The 
P.  0.  box  3707 
Seattle,  wa  98124 

Ola  attn  g.  Keister 

Ola  ATTN  0.  MURRAY 

Ola  attn  G.  hall 
Ola  ATTN  j.  Kenney 

CALIFORNIA  AT  SAN  DIEGO.  UNIV  OF 

P.O.  BOX  6049 
SAN  01  EGO.  CA  92106 

Brcmn  Engikering  Company .  Inc. 
Cummings  Research  Paw 
Huntsville.  Al  35807 

Ola  attn  Romeo  a.  Deliberis 


Charles  Stark  Draper  Laboratory,  Inc. 
S55  Technology  Square 
Cambridge.  Ma  02139 
OIcy  httn  0.  B.  Cox 
OIcy  attn  j.  p.  Gilmore 


General  Electric  Company 
Space  Division 
valley  Forge  Space  Center 
Goddard  Blvd  King  of  Prussia 
P,  0.  Box  8555 


Philadelphia,  Pa  19101  _ 

rvt rv  Attn  M.  h.  SORTNER  jPACE  SCI  lAB 


General  Electric  Company 
?.  0.  Box  1122 
Syracuse.  nY  13201 


General  Electric  Company 
Tempo-Center  for  Advanced  studies 
818  State  Street  (P.O.  Drawer  M 
Santa  Barbara,  Ca  93102 
Ola  attn  uASIAC 
Ola  Attn  don  Chamjler 
Ola  httn  Tom  saprett 
ula  httn  Tim  Stephans 
Ola  httn  warren  S.  xnapp 
Ola  httn  william  mcnamara 
ula  Attn  b.  uambill 

.  1 1/-V  iiTTN  MArv  Qtja/t/tw 


seneral  Electric  Tech  Services  Co..  Inc. 

hmes 

Court  street 
Syracuse,  NY  13201 

Ola  attn  G.  Hillman 


general  research  Corporation 
Santa  Barbara  Division 
P.  0.  Box  6770 
Santa  Barbara.  Ca  93111 
Ola  ATTN  JOHN  lSE  JR 
Ola  ATTN  JOEL  Garbarino 


Geophysical  Institute 
university  of  Alaska 
Fairbanks.  pk  99701 

(All  Class  Attn:  security  officer) 
Ola  attn  T.  n.  Davis  (Uncl  Only) 

01a  attn  Neal  brown  (Uncl  Only) 

Ola  attn  Technical  library 

GTE  Sylvania,  Inc. 

Electronics  Systems  grp-Eastern  div 
77  A  Street 
NEEDHAM.  <V  02194 

OIcy  Attn  marshal  u»oss 


combat  laboratories 
Linthioum  Road  _ 
Clarksburg.  Md  20734 
ula  attn  g.  hyoe 


Illinois,  University  of 
107  Coble  hall 
150  Daveiport  House 
CHAfPAIGN,  IL  61820 

(All  Corres  attn  Dan  .McClelland) 
OIcy  For  K.  Yeh 


Cornell  University 

Department  of  Electrical  Engineering 

Ithaca.  NY  14850 

Ola  HTTN  D.  T.  HARLEY  UR 

O.ECTR0SPACE  SYSTEMS.  INC. 

oox  1359 

RICHARDSON.  TX  75080 

Ola  attn  h.  logston 

Ola  httn  security  (Paul  Phillips) 


INSTITUTE  FOR  DEFENSE  ANALYSES 

jOO  army-navy  Drive 
ARLINGTON.  VA  22202 

ola  attn  j.  m.  aein 
Ola  attn  Ernest  bauer 
Ola  attn  hans  holfhard 

Ola  ATTN  JOEL  BENGSTON 


£51  Inc. 

h95  java  OR ive 

Sunnyvale.  Ca  54086 

Ola  Attn  j.  Roberts 

Ola  HTTN  JAMES  i'VRSHALL 

ula  attn  C.  w.  prettie 
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hSS,  Inc. 
i  Alfred  Circle 
BEDFORD,  HA  01730 

OIcy  Attn  Donald  hansen 


Intl  Tel  &  Telegraph  Corporation 
500  WASHINGTON  AVENUE 
Nutley,  iiJ  07110 

01a  attn  Technical  library 


RAf®  Corporation,  The 
1700  Main  Street 
Santa  Monica,  Ca  90406 
01a  attn  Cullen  Crain 
01a  Attn  £d  Bedrozian 

Riverside  Research  Institute 
SO  west  End  avenue 
new  York,  nY  10023 

01a  attn  vince  Trapani 


JAYCQR 

1401  Camino  Del  Mar 
Del  Mar.  Ca  92014 

01a  Attn  S.  R.  Goldman 


johns  Hopkins  University 
Applied  Physics  Laboratory 
jows  hopkins  Road 
Laurel,  Md  20810 

01a  attn  Document  librarian 
01a  Attn  Thomas  Potemra 
01a  Attn  john  Dassoulas 

Lockheed  Missiles  &  Space  Co  Inc 
P.  0.  Box  504 
Sumiyvale,  Ca  94088 

01a  Attn  Dept  60-12 
01a  attn  D.  R.  Churchill 


LOCKHEED  MISSILES  AM)  SPACE  CO  INC 
3251  hanover  Street 
Palo  Alto,  ca  94304 

01a  attn  martin  Walt  Dept  b2-10 
ola  attn  Richard  G.  johnson  Dept  b2-12 
uIcy  attn  w.  l.  Imhof  Dept  52-12 


Science  applications.  Inc. 

P.  u  BOX  2351 
la  Jolla.  Ca  92038 

01a  ATTN  LEWIS  (1.  LINSON 
Ola  ATTN  DANIEL  A.  HAM.IN 
01a  ATTN  D.  BACHS 
Ola  ATTN  t.  A.  STRAKER 
Ola  Attn  Curtis  \  Smith 
01a  ATTN  JACX  MCDOUGALI. 


Raytheon  Co. 

528  Boston  Post  Road 
Sudbury,  Ma  01776 

01a  Attn  Barbara  Adams 

Science  applications.  Inc. 

1710  Goodridge  Dr. 

McLean,  VA  22102 

Attn:  j.  Cockayne 

Lockheed  Missile  &  Space  Co.,  Inc. 
Huntsville  Research  &  tNGR.  Ctr. 
4800  Bradford  Drive 
Huntsville.  Alabama  35807 
attn:  dale  h.  Davis 


kaman  Sciences  Corp 
P.  0.  BOX  74b3 
COLORADO  Springs.  Co  S0933 
Ola  ATTN  T.  MEAGHER 


III*ABIT  CORP 
10453  ROSELLE 
San  Diego.  Ca  92121 

01a  attn  Irwin  Jacobs 


fl.I.T.  LINCOLN  LABORATORY 
P.  0.  Box  73 
LEXINGTON,  MA  02173 

01a  ATTN  DAVID  M.  TOWLE 
01a  ATTN  P .  WALDRON 
01a  ATTN  L.  LOUGH.IN 

Ola  attn  D.  Clark 


.martin  marietta  Corp 

ORLANDO  DIVISION 

P.  0.  BOX  5837 
ORLANDO,  PL  32805 

01 cy  Attn  R.  meffner 


Physical  dynamics  Inc. 

P.  0.  Box  3027 
Bellevue,  wA  98009 

01a  attn  E.  J.  Fremouw 


Physical  uynamics  Inc. 
P.  o.  BOX  10567 
OAKLAND.  Ca.  94610 
attn:  a.  Thomson 


McDonnell  Douglas  Corporation 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
01a  ATTN  N.  HARRIS 
Ola  ATTN  J.  IMOULE 
Ola  ATTN  OEORGE  iMROZ 
01a  attn  w.  Olson 
01 cY  ATTN  r.  w.  halprin 
OIcy  attn  Technical  library  Services 

mission  Research  Corporation 
735  State  Street 
Santa  Barbara,  Ca  93101 
01a  Attn  P.  Fischer 
01a  Attn  w.  P.  Crevier 
01a  attn  Steven  l.  Gutsche 
01a  attn  D.  Sappenfield 
01a  Attn  R.  oogusch 
OIcy  Attn  r.  hendrick 
01a  attn  Ralph  Kilb 
01a  attn  Dave  Sowle 
OIcy  Attn  F.  Fauen 
OIcy  Attn  m.  Scheibe 
01a  attn  Conrad  l.  Longmire 
01a  Attn  warren  a.  Socueter 

Mitre  Corporation,  The 
P.  0.  Box  208 
Bedford,  ma  01730 

01a  attn  JOHN  Morganstern 
01a  ATTN  G.  HARDING 
01a  Attn  C.  t.  Callahan 


R  4  o  Associates 
?.  0.  box  9695 
Marina  Del  Rey.  Ca  90291 

OIcy  attn  Forrest  Gilmore 
ola  Attn  »yan  Gabbard 
OIcy  attn  william  B.  **ught  jr 
01a  Attn  Robert  F.  lelevier 
01a  attn  william  j.  Karzas 
01a  attn  h.  Ory 
OIcy  Attn  C.  'MacDonald 
Ola  attn  R .  Turco 


mitre  Corp 

westgate  Research  Park 
1820  Dolly  Madison  BLvd 
'MCLEAN,  VA  22101 

01a  attn  w.  hall 
Cla  attn  w.  Foster 

Pacific-Sierra  Research  Corp 
1456  Cloverpield  Elvd. 

SANTA  rbNICA,  Ca  90404 

Ola  attn  £.  c.  field  jr 
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Pennsylvania  state  university 
Ionosphere  research  us 
318  Electrical  Engiicering  east 
University  park,  pa  16802 

(NO  CLASSIFIED  TO  THIS  ADDRESS) 

OIcy  attn  Ionospheric  Research  lab 

Photometrics.  Inc. 

442  Vrrett  Road 
LEXINGTON,  V  02173 

OICY  attn  Irving  L.  kofsky 


reotcLOGY  International  Corp 
75  wiggins  Avenue 
Seoforo,  «a  01730 

OIcy  Attn  W.  P.  Boojist 


TRW  Defense  8  Space  Sys  Group 

One  Space  Park 

keoonoo  Beach,  Ca  90278 

01a  ATTN  R.  K.  PLEBUCH 

01a  Attn  S.  Altschuler 
Ola  attn  D.  Dee 
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